Sub-Poissonian phononic population in a nanoelectromechanical system 
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Population of a phononic mode coupled to a single-electron transistor in the sequential tunneling 
regime is discussed for the experimentally realistic case of intermediate electron-phonon coupling. 
Features like a sub-Poissonian bosonic distribution are found in regimes where electron transport 
drives the oscillator strongly out of equilibrium with only few phonon states selectively populated. 
The electron Fano factor is compared to fluctuations in the phonon distribution, showing that all 
possible combinations of sub- and super-Poissonian character can be realized. 
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Introduction. — Condensed-matter physics and optics 
have often exchanged concepts and ideas, based on the 
common underlying structure of wave phenomena. These 
are essentially based on interference effects which, in the 
case of light, can be derived from classical wave equa- 
tions. However, features exist which cannot be explained 
within a classical treatment, as e.g. squeezed states [l| 
and photon antibunching in resonance fluorescence [2]. 
On the other hand, in mesoscopic physics the focus is 
often on electron transport. The latter is naturally char- 
acterized in terms of current, and much information (e.g. 
carriers charge, process statistics, correlation effects) can 
be extracted from noise Of particular interest are 
then condensed-matter systems in which fermionic and 
bosonic degrees of freedom are coupled, and where elec- 
tron transport induces the emission of nonclassical radia- 
tion. For example, it is known that electronic shot noise 
in a quantum point contact may be source of antibunched 
photons and emission of antibunched phonons from 
a two-level quantum dot is expected when transport is 
characterized by bunching of tunneling electrons [5|. 

Nanodevices where mechanical motion is coupled to 
electric transport constitute in this sense a perfect sub- 
ject for study @. Realizations of these nanoelectrome- 
chanical systems (NEMS) have been obtained e.g. with 
single oscillating molecules 0, semiconductor beams Q 
and suspended carbon nanotubes NEMS are inter- 
esting dynamical systems and are expected to show many 
peculiar transport features ranging from shuttling insta- 
bility to avalanche-like transport (Tlj. Recently, at- 
tention has also been focused on the mechanical proper- 
ties of NEMS, as it appears now that experiments are 
close enough to the quantum limit to test theo- 
retically predicted quantum features in the vibrational 
motion [13[. Of particular interest are those associated 
with the discrete energy states of the oscillator and, in- 
deed, several proposals have been put forward to mea- 
sure discrete number states [13]. Furthermore, it is 
well known that the distribution of oscillation quanta 
(phonons) in NEMS is strongly affected by the trans- 
port of electrons l^, 3], and even the existence of non- 



classical number states induced by tunneling has been 
predicted [17 1. 

In this Letter, we address the behavior of a harmonic 
oscillator coupled to a quantum dot, focusing on the dis- 
tribution of unequilibrated phonons induced by electric 
transport. We show that it is possible to achieve a selec- 
tive population of few phonon states such that the distri- 
bution of the phonon number I displays a sub-Poissonian 
behavior, i.e. Va.1l < (l). At the same time, we con- 
sider the zero-frequency current noise and show that the 
fluctuations of both the phonon distribution and the elec- 
tron current can be enhanced or reduced with respect to 
Poissonian statistics one independently of the other. 

Model & methods. — The system we consider is a gated 
single-electron transistor (SET) coupled to leads and to a 
harmonic oscillator. The SET Hamiltonian is described 
within the standard constant-interaction model for spin- 
less particles. In particular, the charging-energy term is 
He = Ec{n — rig)'^ , where {n — rig) is the effective num- 
ber of electrons on the SET and rig is proportional to the 
charge induced by the gate. The oscillator and coupling 
terms are {h = 1) 



iJph = ujQ b^b + VXujQ {b^ + b){n- Ug), 



(1) 



where W creates vibrational excitations of energy wq- 
The dimensionless parameter a/A defines the strength 
of electromechanical interaction between the position of 
the oscillator and the effective charge on the SET. Such 
a term can be induced by an oscillating gate capaci- 
tively coupled to the dot [8| . The leads are Fermi liquids 
with Hi = ^ ^k,oLc\. ^Cj. ^ and chemical potentials 

Mi,-R ~ Mo ± eV/2, where V is the bias voltage. In the 
limit luq, eV,kBT <C Ec,AE, where AE is the average 
single-particle level spacing, the SET excess occupancy 
is limited to 0, 1 and we can focus on the lowest unoc- 
cupied single-particle level ^. The total Hamiltonian can 
then be written as = en + ffph + Hi + Ht, where 
e = ^ -I- 2Ec{l/2 — Ug) and n = d^d are respectively the 
energy and the occupation number of the single level, 
and Ht — Z)fc,a=L,ii(*ac|,_^c? 4- h.c). Here, ta are the 
tunneling amplitudes, with asymmetry A ~ |iflP/|iLp- 



2 



Being interested in the weak-tunneling limit, we treat 
Ht as a perturbation. It is then convenient to per- 
form a canonical transformation to make the unper- 
turbed Hamiltonian diagonal in the system variables n, /. 
The desired transformation is the Lang-Firsov polaron 
transformation O = UOU^ with U = exp77(& — 6^) and 
r] = \/X(n — Ug) [iBj . The transformed Hamiltonian is 
given hy H = ln-\- ojQb'^h + Hi + Ht, where e = e — Xuiq 
and 



Ht = 



k,a=L,R 



{t^cl^de-^^"'-"^ +h.c.). 



(2) 



In the polaron picture, coherences between states with 
different phonon number can be neglected as far as the 
level broadening 7 induced by tunneling is the smallest 
energy scale into play [3], i.e. 7 <C oja^ksT. In this 
limit, the reduced density matrix p of the SET4-oscillator 
system in the polaron picture is diagonal both in n and 
I, and the dynamics is well described by the rate equa- 



|r'„',.)Pi(0),;' 



,(o)^0(l),ij 



for the 



tions 5tPo(i),i = Ei' i-o{^y 

populations Pn,i = (n,Z|p|n,i/. neic, i^^^^ i^^^a^o) 
are the total rates for tunneling in (out of) the level, and 



\pM. Here, r'/;„^ = r„"' 



I'l 



(3) 



and f{x) is the Fermi func- 
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^ are the 
and V is the den- 



where fa{x) = f{x + e- flo 
tion. The coefficients = 
Franck-Condon (FC) factors 
sity of states of the leads. In the following we assume 
/iQ = ^ — XujQ so that Ug = 1/2 defines on-resonance con- 
ditions. We focus on the regimes of weak (A ^ 1) and 
intermediate (Awl) phonon coupling, where cotunneling 



is negligible out of the Coulomb-blockaded regions [ll | . 

Electronic and phononic expectation values can be 
evaluated in the polaron picture as (O) = Tr[(5p] = 
E„;(n, Z|0|n, l)Pn,i- It is useful to define also a "hybrid" 
average {0)p = Tr[Op]. In terms of (•)p, we can write 

(0 = {l)p + {v')p, (4) 
if) = {f)^ + i{^^l)^+{rf), + {rf)„ (5) 

where we have used the fact that b = b — rj and that p 
is diagonal in the considered weak-tunneling limit. Note 
that for operators like n, which are unchanged by the 
canonical transformation, it is {n)p = (n). 

From the stationary solution of the master equation, 
the phonon Fano factor Fph = Var//(Z) can be directly 
calculated in terms of Eqs.ld]) and The electronic 



Fano factor F = S/2e{I) is evaluated following Ref. [18. 

Results. — We first consider the phonon Fano factor 
Fph. Our main result is that, in the presence of asym- 
metry of the tunneling barriers and specific voltage con- 
ditions, the phonon distribution shows a sub-Poissonian 
behavior Fph < 1 (see Fig. [TJi). 




FIG. 1: a) Phonon Fano factor _Fph as a function of voltage 
V and Tig, at A — 0.1 and A = 0.7. 6) Density plot of Fph 
as a function of A and A for eV = 2ljo, "■g = 0.525 (middle 
of region I in a). Black line: contour _Fph ~ 1. The region 
between the dashed lines encloses the four-state regime (see 
text). In both panels: fcsT — 0.01 a;o, Ec = IOljo- Color 
scales on the top. 



In particular, for ^ < 1 (A > 1) the most favorable 
region of the V — Ug plane for having i^ph < 1 is region I 
(region III) . A sub-Poissonian Fph can also be obtained in 
region II, but only in the limit of very strong asymmetry 
(not shown). For definiteness, in the following we assume 
A < 1 and we focus mainly on region I which, in general, 
is limited by the following conditions: cjq < eV < Sluq 
and 1/2 < Ug < l/2 + ujo/2Ec. 

Here, the phonon Fano factor shows a crossover be- 
tween sub- and supcr-Poissonian behavior as a function 
of A and A (see Fig. As a rule of thumb, a sub- 

Poissonian Fph requires A < 1 and it is favored by strong 
asymmetries. Interestingly, for intermediate values of the 
electron-phonon coupling A « 1 , it is Fph < 1 already for 
asymmetries which are experimentally feasible, A < 1. 

Theoretically, the super- (sub-) Poissonian character 
of Fph is more easily studied in terms of the parameter 
Q ~ VarZ — (Z), being Fph < 1 only if Q < 0. Let us 
write Q = Qp + AQ, where Qp = {P)p - {1)1 - {l)p and 

AQ = 2\{l)p[nl - [2ng - 1)(2 - {n)p)\ (6) 
+ 4A(2ng - 1)(Z(1 - n))-p + \^(2ng - l)2Varn. 

Taking into account that n G {0, 1}, it is easy to show 
that it is always AQ > in region I. Therefore, a sub- 
Poissonian Fph requires necessarily Qp < 0. 

The quantity Qp can be evaluated in terms of the 
phonon distribution in the polaron frame Vi = Fq./ + 
Note that if Vi obeyed a thermal distribution Vi = 
g-/3;a;o(i _ e-/5"o)^ it would be Qp = (e'^"" - 1)-^ > 0, 
while if Vi follows the Poisson distribution it would obvi- 
ously be Qp = 0. On the other hand, when only the two 
lowest vibrational levels are occupied, 'Pq + T'l = 1, it is 

Qp = -{i?p<^- ^ 

A phonon distribution having only the first few states 
occupied is therefore a promising candidate for observing 
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Fph < 1. Indeed, our numerical investigations strongly 
suggest that a sub-Poissonian phonon Fano factor re- 
quires a slender phonon distribution: in fact Fp^ < 1 is 
solely observed at low voltages, in the presence of asym- 
metry and preferably off-resonance, which are all condi- 
tions which favor a narrow Vi ■ 

Remarkably, for A = 1 the exact stationary solution 
of the rate equation satisfies the condition Vq + Vi — 1 
always in region I. This fact can be understood observing 
that in region I the only energetically allowed transitions 
which increase the phonon number arc (0, 1) — > (1, / -I- 1), 
sec Eq.Q. As a consequence, at low temperatures 
kgT <C Wo, excited phonons states can only be popu- 
lated via a series of subsequent tunneling events such as 
(0,0) -> (1,1) -> (0,1) -> (1,2). However, for A 1 
the transition (1, 1) (0, 1) is forbidden because the FC 
factor X^^ = e^^{l — A)^ vanishes. In this case, the dy- 
namics of the system is frozen to the states with I < 1 and 
it can be solved exactly by considering the reduced four- 
state model represented in Fig. [2li. Within this model, 
we obtain 



Ph|A=l 



1 



3 — 8n„ 



2A{3-8ng + 5nl) 



i2 + A)[n^{A + 2)-4ng + i] 



which gives e.g. Fph = 0.96 for Ug = 0.54 and A = 0.1. 

The four-state model of Fig. [2}i is often a good ap- 
proximation of the full numerical solution also for A 7^ 1 . 
Within this model, we derive an analytical expression 
for Fpii which, however, is too long to be reported here. 
This approximates the exact numerical result with an er- 
ror smaller than 1% in all region enclosed between the 
dashed lines in Fig. and therefore it allows us to in- 
vestigate analytically the crossover between super- and 
sub-Poissonian phonon Fano factor. In particular, wc 
find that it can be Fph < 1 only for values of the electron- 
phonon coupling smaller than a certain critical value Acr 
which, up to order [ug — 1/2)^, is given by 



Ac 



4(ng - 1)2 + 2A(n2 - 4ng -t- 2) ■ 



(7) 



This equation describes very accurately the upper part 
of the contour line i^ph = 1 in the phase diagram Fig. [2), 
deviating from the exact result only for A 1. 

The existence of a maximum critical value Acr < 1 
can be understood qualitatively considering the limit of 
strong asymmetry A 1. In this case, it is (n) « 1 
and Varn sa 0, and from Eq.® one obtains directly Q 
(Z)p[2A(ng — 1)2 — where we have used Qp = —{l)^p in 
the four-state model. A strong electron-phonon coupling 
is thus unfavorable for Q < in two respects; on one 
hand, it increases the weight of the positive term oc (ug — 
1)2; on the other one, it is well known that Vi=o — * 1 as 



A is increased [la], so that (l)p 0. 



We can conclude that the sub-Poissonian Fano factor 
is induced by a phonon distribution Vi in which only 
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FIG. 2: a) Set of states included in the four-state model. 
States are labeled as (n, /), arrows represent the relevant tran- 
sitions in region I for ksT ujq. b) Vi (black) vs. Vi (gray) 
for eV = 2ujo, Ug = 0.525, fcsT = O.OIljo, = Wujo and 
A = 0.9, A = 0.01. 



the first few phonon states are populated and yet the 
occupation probability of the excited states is comparable 
with the one of the ground state. We refer to such a 
situation as a selective population of the phonon states. 

Finally, we remind that Vi is the phonon distribution 
in the polaron picture. The intrinsic phonon distribution 
is Vi = J2n Pri,i where Pn,i = (n, l\p\n, I) and p = WpU 
is the density matrix in the original picture. Note that, 
in terms of Vi, the average phonon number reads (l) = 
'^ilVi, and similarly {P) = "^il^Vi- We can evaluate 
Vi taking into account that Pn^i = ^n,i' , where 

-^n — \ {TT'J\U\n,l')\'^ are generalized FC-factors. Such a 
relationship is a consequence of p being diagonal in the 
weak-tunneling limit. Comparing Vi and Vi, it is clear 
that one can speak of selective population only in the 
polaron picture (sec Fig.[2t). However, what is important 
is that, in the presence of a selective population of Vi, the 
intrinsic phonon distribution Vi shows a sub-Poissonian 
behavior, signaled by i^ph < 1- 

Up to now, we have considered solely the character- 
istics of the phonon distribution induced by tunneling. 
However, it is well known that the transport properties of 
the system are in turn strongly affected by phonons. Sig- 
natures of this interplay are especially visible in the cur- 
rent Fano factor, which is very sensitive to the electron- 
phonon interaction [ill . [l9| . For example, a giant en- 
hancement of F has been predicted as fingerprint of 
strong electron-phonon coupling [llj. Here, we consider 
intermediate coupling and wc focus on the study of the 
(sub-) super-Poissonian character of F with respect to 
the one of _Fph. Interestingly, in region I all the possible 
combinations of F, i^ph ^ 1 can be obtained by tuning the 
asymmetry A and A (see Fig. [5]i). This is possible be- 
cause the super- and sub-Poissonian character of F and 
-Fph have different physical origins. In fact, while i^ph < 1 
presumes a selective population of phonon states, F > 1 
is induced by a bunching of tunneling events [20| . 

A simple explanation of this mechanism can be given 
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FIG. 3: a) Phase diagram of the possible combinations of 
f ph in region I [eV = 2ujo,ng = 0.525) depending on A 
and A (see text for discussion). h) Plots of Fph and F 
as a function of V for Ug = 0.525 and A = 0.6, A = 0.01 
corresponding to the asterisk in panel a for eV = 2uJo. Boxes: 
exact numerical solutions; lines: four-state approximation. 
In both panels fcsT = 0.01 o^o, Ec — lOu^o. 



in terms of the four-state model of Fig. For A 
the electronic Fano factor can be written as: 



1, 



4 + ^2 4: + 2A + UA'^ + 9A^ - 



(2 + Ay 



{1 + Af{2 + Af 



For A = 1 the system behaves as a spin degenerate sin- 
gle level so that it is always F < 1 [2lj. For Awl, 
the transitions (0,0) ^ (1,1) and (0,0) ^ (1,0) act as 
two competing transport channels, whose relative weight 
is determined by the ratio X^^/X'^^ = A. It follows 
that for A < 1 the state (1,1) is a trap state and blocks 
the transport through the other more conducting chan- 
nel (0,0) <-> (1,0). In the presence of asymmetry, such 
a dynamical channel blockade [lOl leads to bunching of 
tunneling events and to super-Poissonian current noise 
F > 1. However, as the difference between the two com- 
peting transport channels is fairly weak, F is only slightly 
above 1 (see Fig. |3f)). The same mechanism occurs for 
A > 1 but, in this case, it is transport through the excited 
state that is blocked by the occupation of (1,0). 

It is then clear why super-Poissonian current noise and 
sub-Poissonian phonon distribution can occur simultane- 
ously only for A < 1 when the trapping mechanism re- 
sponsible for F > 1 also favors the selective population 
of the phonon states. Viceversa, for A > 1 the occupa- 
tion of the vibrational ground state is strongly favored 
since (1, 0) is the trap state, and the phonon distribution 
is mainly super-Poissonian. Note that outside region I 
the four-state model differs considerably from the exact 
results, which exhibit F < 1 and i^ph > 1 as expected 
from a fermionic and bosonic system, respectively. 

Finally, a comment is in order. A suppressed phonon 
Fano factor i^ph < 1 has recently been predicted for an 
oscillator driven by a superconducting SET in the limit 
7 Wo, and this has been interpreted as signature of 
a number-squeezed state [17| . In our case, instead, we 
obtain a sub-Poissonian distribution without squeezing. 



This is ultimately a consequence of the loss of phase in- 
formation in the weak-tunneling limit 7 <C |22| . 

In conclusion, we have shown that a sub-Poissonian 
phonon distribution can be achieved in a nano electrome- 
chanical system when tunneling induces a selective pop- 
ulation of few phonon states. In addition, we have con- 
sidered the electronic noise and we have found different 
combinations of sub- and super-Poissonian electron and 
phonon Fano factors depending on the asymmetry and 
on the strength of the electron-phonon coupling, 
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